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Abstract. We deal with an initial-boundary value problem for the gen¬ 
eralized time-dependent Schrddinger equation with variable coefficients 
in an unbounded n-dimensional parallelepiped (n > 1). To solve it, the 
Crank-Nicolson in time and the polylinear finite element in space method 
with the discrete transparent boundary conditions is considered. We 
present its stability properties and derive new error estimates 0{t^-\-\ h\^) 
uniformly in time in space norm, for n > 1, and mesh space norm, 
for 1 < n < 3 (a superconvergence result), under the Sobolev-type as¬ 
sumptions on the initial function. Such estimates are proved for methods 
with the discrete TBCs for the first time. 
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1 Introduction 

The linear time-dependent Schrddinger equation is the key one in many phys¬ 
ical fields. It should be often solved in unbounded space domains. A number 
of approaches were developed to deal with such problems using approximate 
transparent boundary conditions (TBCs) at the artificial boundaries. 

Among the best methods of such kind are those using the so-called discrete 
TBCs remarkable by the clear mathematical background and the corresponding 
rigorous stability results in theory as well as the complete absence of spurious 
reflections in practice. They first were constructed and studied for the stan¬ 
dard Crank-Nicolson in time finite-difference schemes, see m and also |2I3) . in 
the cases of the infinite or semi-infinite axis and strip. Later families of finite- 
difference schemes with general space averages were treated in mm- In par¬ 
ticular, they include the linear and bilinear FEMs in space. 

In this paper, we consider the Crank-Nicolson-polylinear FEM in an un¬ 
bounded n-dimensional parallelepiped (n > 1), present results on its stability 




with respect to the initial data and a free term as well as on exploiting the dis¬ 
crete TBCs and mainly derive the corresponding new error estimates 0 (t^ + 
uniform in time and in norm (for n > 1) and mesh norm (for 1 < n < 3) 
in space under the Sobolev-type assumptions on the initial function. The latter 
estimate is a superconvergence result. Such estimates are proved for the meth¬ 
ods with the discrete TBCs for the first time. Importantly, the error estimates 
contain no mesh steps in negative powers like for other approximate TBCs, see 
m, that is one more advantage of using the discrete TBCs. 

2 The IBVP and numerical methods to solve it 

We deal with the initial-boundary value problem (IBVP) for the time-dependent 
generalized Schrodinger equation with n > 1 space variables 

ihpDtil) = T-Lip dw{BVil))+ Vtl) on 7Tx(0,T), (1) 

Man = V'|t=o = on 7T. (2) 

Hereafter ip = ip{x^ t) is the complex-valued unknown wave function, i is the 
imaginary unit and > 0 is a physical constant. The x = {xi ,..., Xn)-depending 
coefficients p,V G L°°{II) and the nx n matrix B G L°°{n) are real-valued and 
satisfy p{x) > p > 0 and B{x) > B/ > 0 on 7T, where I is the unit matrix 
(whereas V can have any sign in general). Here U := R for n = 1 or, for n > 2, 
7T := R X TTy is the infinite parallelepiped, with TTy := (0, X 2 ) x • • • x (0, Xn)- 
Also Dt = ^ and Di = ^ are the partial derivatives, and the operators div 
and V are taken with respect to space variables. 

We also assume that, for some (sufficiently large) Xq > 0, 

p{x) = Poo, H(x) = diag(Bioo,---,S„oo), V{x) = Voo for |xi| > Aq, (3) 

where diag {Bioo, ■ ■ ■, Bnoo) is the diagonal matrix with the listed positive diag¬ 
onal entries. More generally, it could be easily assumed that p, B and V have 
different constant values for xi < —Xq and for xi > Xq. Let Xi > Xq, and 
f2 = fix = (—Ai, Ai) for n = 1 or 17 = fix = (—Ai, Ai) x TTy for n>2. 

We consider the weak solution ip G C'([0, T]; i7g (7T)) having Dtip G C([0,T]; 
L^[n)) and satisfying the integral identity 

ih{Dtip{-,t),(p)L 2 ,,^jj-^ = £„{iP{-,t),(p) for any <p G H^{n), on [0,r], (4) 

and the initial condition ip\t=o = ip^ G Hl{n). Hereafter we use the stan¬ 
dard complex Lebesgue and Sobolev spaces (and subspaces), the weighted com¬ 
plex Lebesgue space L^’^(G') endowed by the inner product (w, ip)l 2 ,p(( 3 ) := 
(pw, (p)i 2 (( 5 ) and the "H-related Hermitian-symmetric sesquilinear form 

£g(w,(p) := ^(BVw, V<p)i 2 (G)-f (Hw, ■p)l 2 (g), with G = 7T, 17, etc. 

We define a non-uniform mesh in xi on R containing the points ±Ai and be¬ 
ing uniform with a step 0 <hi< Xi outside [—Xi+hi,Xi — hi] D [—Ao,Ao]. We 


also define non-uniform meshes in X 2 , ■ ■ ■ ,Xn respectively on [0, X 2 ], ■ • ■, [0,X„] 
(containing the ends of the segments). They induce the partition of 77 into finite 
elements that are rectangular parallelepipeds without common internal points. 
Let |li| be their maximal diagonal length. Let Sh{n) be the (infinite-dimensional) 
subspace of functions in Hl{n) that are polylinear over each element. Clearly 
Sh{n) C (7(77) n L'^{n). Let Sh be the restriction of Sh{n) to fl. 

Let ujj^ be the non-uniform mesh 0 = to < ... < Im = T with steps := 
tm-tm-i- We put Tmax := niaxi<m<M Tm and fm ■= for 1 < m < M-1 

and To := We define the time mesh operators 

_ ■y^m ■y^m—1 y^m+1 y/'m ■ym —1 i ym 

dtY^:= -, dtY^:= -^-, StY^ := -^-. 

We introduce the Crank-Nicolson-polylinear FEM approximate solution F: 
Um —>■ Sh{n) satisfying the integral identity 

F)L^’P{n) — f) for any tp g Sh{n) and I < m < M, (5) 

compare with (j4]), and the initial condition '7'|t=o = F^ & Sh{n), where is an 
approximation for 

Let be a conjugate linear functional on Sh{n) that we add to the right- 

hand side of ([5]) to study stability in more detail and to derive error estimates. 


Proposition 1. Let P^{^) = {F"^, F)L'^(n) with F™ g L'^{n) for 1 < m < 717. 
Then there exists a unique approximate solution F and the following first stability 
bound holds 


M 




( 6 ) 


We introduce also the “energy” norm such that 

\\w\\H+vp-n ■= i^n{w,w) +v\\w\\l2,f,^^n) > ^^\\w\\l^,p(n) for any w g Hl{n), (7) 

with some real numbers v and 5 > 0. Inequality © is knowingly valid for v so 
large that \B_Xq -|- F(a;) -|- (7 — S)p(x) > 0 on 17 with Aq := J2k=2(^)^ (here 
Ao = 0 for n = 1). We define also the corresponding dual mesh depending norm 

ll^llm-bm := ™ax \{w,ip)n\<c\\w\\H-^n), 

^ > <peSh(7I):||¥>||K+«p;77 = l 

where {w, (p)n is the conjugate duality relation on H~^{n)x 77q (77) and 
[77g(77)]*. Hereafter c and ci are generic constants independent of the meshes, 
any functions and T whereas cq denotes absolute constants (fixed numbers). 






Proposition 2. Let CM = {F"",V>)n with F™ € H-^{n) for 1 < m < M 
and F° G H~^{n) be arbitrary. Then there exists a unique approximate solution 
T and the following second stability bound holds 

o<m<M" "H+vp,n - \\ wn+vpin 

M |.| 

+4 (!^ Ilf-IU-.m + + 4 ||F<>||;;‘> ^ (8) 

m—1 

Method (O cannot be directly used in practice because of the infinite number 
of unknowns at each time level. Nevertheless it is possible to restrict the method 
to 17 provided that G Soh ■= G Su] g}{x) = 0 on 17 \ l7o}, where f7o := 
L2xi-hi,X2,...,XnJ ZJfj is uniform with the step t = ^. Let both assumptions 
be valid up to the end of the section. 

By definition, the discrete TBCs are conditions at the artificial boundaries 
xi = ±7fi allowing to accomplish the restriction (they are non-local uiX 2 , ■ ■ ■ ,Xn 
and t). To write down them explicitly, for clarity, we confine ourselves by the 
case of the uniform mesh in Xk with the step hk = for 2 < k < n, and define 
the related well-known direct and inverse discrete sine Fourier transforms 

^ Pj sin l<q<Jk- 1, 

Jk , Ok 

J = 1 

Pj = {Pk^P^'^), ■■= E P^'^ 1 < J < Jfc - 1. 

9=1 

The related eigenvalues of the ID linear in Xk FEM counterparts of the operators 
—D^ and the unit one (for zero Dirichlet boundary values at Xk = 0,7ffe) are 



Denote by the internal part of the introduced uniform mesh in i7y and define 
the related mesh inner product 

.72-1 J„-i 

(C7,1TC^ := E E for n>2 

92 = 1 = l 

or set (D, := UW* for n = 1, where W* is the complex conjugate for W. 

Recall that the discrete convolution of mesh functions It, Q: —>■ C is given 

by {R * Q)™ := E™=o for 0 < TO < M. 

Proposition 3. The restriction 'I\q of the above approximate solution obeys 
the integral identity on 12 

ih(dtT^,g^)L2.Pin) = Cn{-stT'^,T) 


( 9 ) 




for any ip G Sh and 1 < m < M, and the initial condition 


flt^o = e Sh. 

Here ’ • ■ • ’ ^™Li-±Xi} ® vector-function. 

The operator iSref in the discrete TBC has the form 


( 10 ) 




-1 


ai? ... 




( 11 ) 


for any <T: ujp^ x ujff —>• C such that = 0, where €>"* := ..., 

:= and q = (g 2 , ■ • ■, Qn)- Here the kernel i?q 


be 


computed by the recurrent formulas 


yiryi _ -< 'ryi _ > 

i?q = Clq, R\ = -ClqKqPq, R'T = -Kq/Xqi?™-m > 2 , 


'-iq? '-iq'^qr^q^ -^q 

with the coefficients defined by 




Clq — 




(arg aq)/2 


5 '^q 


= _e--g“., = (12) 


iq = 2 aq + i/iiflq ^ 0 , arg ttq S ( 0 , 27r) , / 3 q = 2 Reaq + i/ii|aqp, 


Kx 


On — 


Rioo/i2 2B 


Bioo 


1,(2) !,(»), 
^O, , , Ag„ 


3 

2/9o 


^92 

The next lemma is important to prove stability results for method (jH]), (flUll . 
Lemma 1. The operator 5™£ satisfies the inequalities IMS 


^ > 0, Im^ {ihdt + vstW)^ _ r > 0 


Z=1 


/ = 1 


on ujJj, for any <P: cvp^ x co]^ —>■ C such that = 0 and v > — ^ (see (|31 ) )■ 

Let £'^{p) be a conjugate linear functional on Sh that we add to the right- 
hand side of © to study stability in more detail. 

Proposition 4. Let e^{p) = with F™ e L^{n) forl<m<M. 

Then the solution to ®, (nni) is unique and satisfies the first stability bound 


max 

0<m<M 


2 ^ 
m—1 


IIF" 


Il2,i/p(£2) 


T. 


(13) 


We introduce the “energy” norm on 17 such that 

\\w\\n+i,p.a ■= ^oiw,w) -G V \\^w\\l2(^si) > 


( 14 ) 











foranyw S := w|(_Xi,Xi)xa 77 j = 0} except for w = 0,and some 

real number v > —In particular, for v so large that ^^Xq+V{x)-\-vp{x) > 0 
on 17, m is valid. Define also the respective dual mesh depending norm 

^ ' ¥>6Sh: ||vllw+®p; « = 1 


where {w, p)n is the conjugate duality relation on H ^(17) x 77^(17). 
Proposition 5. Let n{ip) = {F'^,‘p)n with F™ e 77-i(f7) for 1 < m < M 
and S FI~^{f2) be arbitrary. Then the solution to ([5]), (1101) is unique and 
satisfies the second stability bound 


max 

0<m<M 




'H-\-vp; Q 


< 


\p^ 


'H+iip-, Q 


M I „ I _ 

+4 y; (M + ||af”||„_.,„,)r+4 ||F“ 

m—l 


H-\Q) ■ 


(15) 


Propositions [T]-[5] and Lemma [T] in the quite similar cases of the semi-infinite 
77 = (0,oo) (n = 1) and 77 = (0, oo) x (0,X2) (n = 2) were proved respectively 
in [3] and mm (see also m). where families of finite-difference schemes with 
space averages depending on a parameter 0 were treated covering, in particular, 
the linear and bilinear FEMs (for 0 = i). For the presented improvement in 
formulas (HU, see also [14]. For n = 1, the results are as well particular cases 
of those from [13] (given specihcally for general FEM). The case n > 3 can be 
treated in the same manner as n = 2 (for such an exampl e, see [Til l. 

The numerical results for the method can be found in 


3 Error estimates 

Let condition 0 be valid and aw be the elliptic projection of w € Ft^ill) onto 
Sh{n) such that 

Cn{(Jw,(p) +v{aw,ip)L2,p(^n) = Cn{w,(p) + v{w,(p)L2,p^n) (16) 

for any ip G Sh{n). Note that aw exists and is unique. We also assume below 
that B G W^’°^(II) and then the following error estimate holds 

||ui - awllL 2 .P(jj) < c\h\^\\{'Hp + 7)w||l2.p( 77) for any w G H‘^(n) D Hl{n). (17) 

We consider 77p := ^77 as an unbounded operator in L‘^’P{n) with VifHp) = 
F[‘^{n) (~l F[l{n). Assume below that G 
Proposition 6. The following first error estimate holds 

-|-c(l -I- T) {r^axll’^pV’°llL2,p(77) + |liP(l|77p'0°||L2,p(77) -|- ||V’°IIl2,p(/7)) } ■ (18) 

Here afP can be replaced by . 






Proof. 1. For y G L^{0,T), define the average (the projection on the time mesh) 
[y]^ := ^ y{t)dt, 1 < m < M, and notice that 

M 

IIM^IIbt™ < / \\u{-,t)\\Bdt, (19) 

m=l "'0 

where || • ||b = || • Wl^.p^h) or || • \\-H+vp-,n and u S L^(0,T;B), and 

\[yr-Sty"^\<coTi[\DM^, l<m<M, for yelF^’i(0,r). (20) 

2. Applying [•] to identity (0]), we get 

ih{[Dtilj]"',‘f)L^,p(^n) = for any e i?Q(7T) and 1 < m < M. ( 21 ) 

Then for any rj: ujfj —>• Sh{n) from identities ([5]) and (1^ it follows that 

ih(dt{'P - ?7)™, yf)L'^.P(n) - - ??)"", p) 

= imiDtip] - dty)"',(p)L^.P(n) - ^n{{[ip] - stv)"',ip) for any ip G Sh{n). 
Let 77 ™ := cr'0™, 0 < m < M. By identity (fTBl) and [Dty] = dty we get 
ih{pdt{d' - p)L^n) - Cnisti'P - ctt/i)™, (/?) 

= *ft(p([AV']™ - [Ao’V']™), 

+D(pst ('0 - p)L^n) =■ (F, p)L^{n) for any 1 < m < M, ( 22 ) 

where (after rearranging the summands) 

F = -{[Htp] - stHtp) + ihp[Dt{i> - atp)] + vpst{ip - atp). (23) 
Let now on = 1. Proposition [I] together with (1191) lead to the bound 

M 

njnax ||(!F-crV')’”llL2.P(77) < \\F° - a'llj°\\L2.f,(^n) +‘^Y1 II^™IIl2,i/p(77) 

0<m<M ^ ’ 

M 

< ||<F° - cr?/;°||L 2 ,p( 77 ) + 2 ^ II {[Hptp] - StT-Lpip) ||L2.p(77)'rm 

m—1 

+ 2 / ||A(V'-<r^)llL2.p(77)C^i + 2|'0|r max II (■!/)-o-!/;)”" 11^2,p(77). 

The formula ip — F = ip — aip — {F — aip) and estimates (fTTll and (120)1 imply 
max \\{iP - F)^\\L2,p(^n) < \\F° - aip^WL^.P^n) 

+c{r-Lx^ \\D^npip\\L2.P(n)dt+\h\^ WDtCHp + v)ip\\L2.P(n) dt 

+(1 + |f;|r)|/i|2 max IK-Hp + f))V'™||L 2 .p(^)|. 

0<m<M ^ ' J 


( 24 ) 


Under the above assumptions, the solution to problem ([T]), ([2]) satisfies the 
bound 

^ ll(A^V')|t=o||L2,p(77) = ||Hpl/’°||L2,p(77), 0 < fc < 3, (25) 
and the property D^-ip = {—iT-Lpfi!) for 1 < 1 < fc. Therefore 

0^ax^ll(V'-^')’"l|L2,P(7T) < ||'f°-Cr?/’°||L2,p(77) + C {Tr^axll'^pV'° l|L2,p(i7) 

+r|hn|H>°|U2,p(;7) + (1 + \v\T)\h?{\W°\\L^,,in) + |fi|||V'°llL^..(77))} . 

Note that for '0 = 0 the estimate is simplified. 

The following multiplicative inequality holds 

ll'^pV'°lli 2 ,p( 77 ) < ||Hp+V°IIl 2 ,p( 77 )||Hp“V°I|l 2 ,p(/ 7 ) for ^ = 1 , 2 . (26) 

Using (II3 for w = ip^ together with (1^51) for Z = 1 , we complete the proof. □ 

Corollary 1. Let ip°{x) = 0 for |a:i| > Xq, S Soh andufj be uniform. Then 
for the solution to the following first error estimate holds 

n max \\{ip - 
0<m<M 

+c(l +T) {T^||Hp'!/'°|lL 2 ,p(f 2 ) + \h\'^{\\'Hl'tp°\\L 2 .P(n) + |IV'°IIl 2 ,p(j 7 ))} . 

The result immediately follows from Proposition | 6 l Notice also that, for 1 < 
n < 3, for the interpolant in Squ for the following error estimate holds 

||s^° - ' 0 °||L 2 ,p(f 2 ) < c\hf {\\'Hp 1 p°\\L‘^.p(^n'j + ||■^/'°||L 2 ,p(J 7 )), 

thus one can set simply 1 ?'° := sijj^. Other possible choices of S'® with the same 
error estimate, for any n > 1 , are the L'^^Hq) (possibly with a weight like p) 
projection of onto Squ or its elliptic projection onto Squ like (fTBl) (with 7T 
replaced by l7o and any tp € Soh)- 

Remark 1. Importantly, all the above results can be essentially generalized rather 
easily. For example, in the 2D case, the problem in an unbounded domain 7T of 
general shape with smooth boundary and several half-strip-like outlets to infin¬ 
ity can be treated by using combined linear triangle elements inside 77 except 
outlets and bilinear rectangular elements inside outlets. 

Let B e W^'°°{n) and p, U G W^’°°{n). 

Proposition 7. Let Tlptp^ G HQ{n), 1 < n < 3 and := sfP. Then the 
following second error estimate holds 

n^iax ||si/'™ -lf™||w+f)p;77 < c(l -f T)|r^^,,(||77^V’°l|w+t)p;77 + WUlif^Wn+vp-.n) 

0<m<M 1. 

+ |hP(||77pV’°|lL2,p(77) -I- ||'0°||j^2,p(77))|. 

Here sip'^ is the interpolant in Sh{n) for 0 < m < M. 


( 27 ) 


Proof. 1. Let first n > 1. Inequality o implies 


\\pw\\H-\n)^^ lk||L2,p(^) < (5 ^^^\\w\\n+vp-,n- (28) 

Then setting c := 1 + we also get 

11'^’*^ I Iff” ^(77) — IK'^ + '^^)'^llff”^(ff) + l'^ll|p^llff”i(ff) ^ c||ui||-H+i}p; ff • (29) 

For y G L^(0,T), define two more averages (projections on the time mesh) 

[y]T ■= [ yit)dt, (y)"" ■= ^ f y{t)em{t)dt, 1 < m < M - 1, 

2'^m Jtm-l 

[y ]2 ■= [y]\ (y)° ■= — [ yit)eo{t)dt, 

■’’l Jo 

where emit) is the “hat” function linear on all segments and such that 

Cm (tm ) = 1 and Bmiti) = 0 for all I ^ m. 

Notice that the following relations hold 

4br = (A2/)™, (ftsty)™ = [A?/]r, l<m<M-l, (30) 
-y(0) = y(A 2 /)°, (sty)^-2/(0) = y[Ay]\ for y G VF^’^(0,r), (31) 

M—1 .T M-1 „T 

ElK«ni^™</ \H-Mdt, EIIM™I|t,«</ \\ui;t)\\dt, (32) 

m=0 m=0 do 

where || • || = || • \\L^.p{n) and u G L^’^(i7 x (0,T)): and 

liy)""-[y]T\<coTm[\Dty\]^, l<m<M-l, for y G IF2’^(0, T)- (33) 

2. Let first G Shin) be arbitrary. We go back to the error identity ([2^ . 
Applying Proposition [5J now we get 

ninax ||(tF-cr?/>)™||«+fip;77 < - af}°\\-H+vp-,n 

0<m<M 

M-1 M 

+4 ^ A + 4|D| ^ + 4||F°||^-.(^), (34) 

m—0 m—1 

where the right-hand side is slightly transformed and F is given by (l23)) . We 
introduce the decomposition 

F = Fr + pFh, Fr :=-il'H'ip]-St'H'tp), Fh ■■= i[Dtitj} - cri/;)] + vstitp - afj) 
as well as set := 0 and = iDtif) — a'if)\t=o + ■ 0 ( 2 /'° — crijj^). 


Applying sequentially relations (1^ . (1501) . (1^ . (1551) . (fTOl) and (15^ . we obtain 


M 


M-1 




< 


Sr := Ir)| ^ \\F:r\\H-\n)rn. + E 

m—1 m—0 

M M-1 

c(l^l E \\i{'^]-St'tp)'^\\'H+vp-,nTm + E “ [A^/’]2)™||w+«p;7I'rm 


+ ll(M-stV') \\n+vp-,n^ 

M M-1 

< CCor^ax(l^^l E \\{\Dt^\r\\H+vp-nT,n+ E II [l-^tVl]™ ||w+Cp; 


m—1 


m—1 


+ ll[l-D?V’l]^l|w+fip;/7 


< 


CCoT-max{ll(A^V’)|t=o||w+*p;77 + J (|r) 11| £>? V-llw+fip; i7 + 21| £>*^11 W+fip; h) dt}- 

(35) 

The left inequality (pS)) implies 


M 


M-1 


< 


:= |D| E \\pFh\\H-\n)rm + WpF^Wn-^n) + E IIM^ril/^-(77)^™ 

m=l m—0 

M M-1 

ri/2(|t)| E iii"riu=-(77)r™ + ii£^°iil2..(^) + E ii^^niL^-mr^), (se) 


m—0 


and further the error estimate m leads to 

I|£'°I|l2,p(7T) < II A(V' - 0-?/')|t=o||L2.p(77) + |^|||'0° - 
< c|/ip(||A('Hp + '0)'i/'|i=o||L2.p(77) + I^IIK'Hp + 'y)V'°||L 2 ,p(^)). (37) 

Applying sequentially relations (l30l) . (1^ . (ITTl) and (1321) . we also obtain 

M-1 

E \\^tF!l'\\L^.P(n)Tm 

m—0 

M-1 

< E (IK^t “ CrV’))™llL2.P(77) + l^ll|[A(V' - 0-V')]™llL2,P(77))Tm 

m—0 

M-1 

< c\h\^ E {WiFU'Fp + v)ilj)"'\\L^.p(n) + l^|||[A('Hp + D)V’]rilL=''>(i7))Tm 

m—0 

<c\hf ( (||T)i ("Hp +'0 )^||l2,p(7j) + InlllA('^p + f^)'*A||L2.p(77)) (38) 

^0 


Inserting into (IMl) all the estimates dsni)-® together with the estimate for 
Sm=r \\Ph^\\L'^'i’(n)Tm used in the preceding proof in (l24l) . we derive 

max IKiIr - (Tij))'^\\-H+vp-, n < - cr'tp^Wn+vp^n + 45^ +45/1 

0<m<M 

< cTLx{llP?V')lt=o||w+«p;77 {\v\\\D^tlj\\n+vp-n + \\D^'tp\\n+vp-,n)dt'^ 

+c\h\'^Uv\ f \\Dt{'Hp + v)'ilj\\L^,p(n)dt + v^T max ||('Hp + ■O)'0 ™||l2,p(77) 

'■ JQ 0<t<T 

+ ||i4i('Hp + v)llj\t=o\\L^.p{n) + I^IIK’^p + ^')V'°IIl2,p(77) 

+ J {llDtiy-p + v)'>P\\L^,p{n) + \'v\\\Dt{'Hp+v)'>p\\L2,p(n)) dtY 

Once again for '0 = 0 the estimate is essentially simplified. 

Under all the above assumptions, the solution to problem (P), (P satisfies 
the bound 

0 ™^ llA^V’ll'H+i)p;77 < \\{Dtip)\t=o\\'H+vp-,n = W'Hp'tp^Wn+vp-, n, 0 < fc < 3. 

This bound and (051) and the property D^tp = D^~\—i'Hpy'ip, 1 < I < k, imply 
max |l(if - crr/>)'"||w+fip;/7 < \W° - fTip°\\n+vp; n 

0<m<M 

< c(l + T)^T^^^{\\'Hl'(p°\\-H+vp-, n + \\'Hl4’°\\n+vp-, n) + \h\‘^ ^ II'^pV’°IIl2.'>(J7) |- 

k^O 

(39) 

3. Let now 1 < n < 3 and = s'0°. Similarly to [9], Lemma 5.1 for n = 1 
and m. Theorem 2.1 for n = 2 and 3 (see also 0), the following elliptic FEM 
error estimate holds 

n 

llsui - (Jw\\-H+vp-n < c|/ip (S \\DpDqw\\L^(^n) + X! 

P#9 P=1 

< ci\h\'^(WHpwWn+vp- n + \\w\\n+vp-, n) 

for any w G 'D{%p) such that Hpiv G Hl[n), taking into account the above 
regularity assumptions on B,p and V, where the first sum is taken over all p 
and q from 1 to n excluding p = q and disappears for n = 1. This estimate allows 
to pass from to the final estimate (1771) by the triangle inequality together 
with inequalities (Ull) and 

ll’^pV'°llw+fip;i7 < IK’^P + ^)’^p'0°I|l2,p(77)||Hp'0°||l2,p(77), Z = 0, 1. □ 

Corollary 2. Let %l/^{x) = 0 for |a;i| > Xq, 1 < n < 3 and = sip^ on 17 and 
Lufj be uniform. Then for the solution to the following second error 


estimate holds 


max IlsV''" - 'I'n\n+vp-n< c{l+T){TW\-Hy\\n+vp ;0 + \\ny\\u+vp-n) 

+ I^P(ll’^p^°IU2.p(n) + ||'0°||L2,p(f2))|. 

The result immediately follows from Proposition [71 

Note that the norm || • \\u+vp-,n is equivalent to || • \\h^(o) and ||s • ||_f/i(j 7 ) is 
actually the mesh counterpart of the latter norm. 
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